function sse = FitGauss_lin(p, xdata,ydata )
        FittedCurve = gauss_lin(p,xdata);      
        ErrorVector = FittedCurve' - ydata;
        sse = sum(ErrorVector .^ 2);
end